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Abstract. We present recent results on generalized parton distributions from 
dynamical lattice QCD calculations. Our set of twelve different combinations 
of couplings and quark masses allows for a preliminary study of the pion mass 
dependence of the transverse nucleon structure. 

1 Introduction 

The investigation of the nucleon structure in the framework of QCD is a major 
task in today's particle physics. Generalized parton distributions (GPDs) have 
opened new ways of studying the complex interplay of longitudinal momentum 
and transverse coordinate space as well as spin and orbital angular momentum 
degrees of freedom in the nucleon. The twist- 2 GPDs H, E, H and E of quarks in 
the nucleon are defined by the following parameterization of off-forward nucleon 
matrix elements 
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Figure 1. Graphical representation of GPDs as part of a scattering amplitude 

for the helicity dependent axial vector case pQESIH]. Wilson lines ensuring gauge 
invariance of the bilocal operators are denoted by U. Here and in the following 
we do not show explicitly the dependence of the GPDs on the resolution scale 
Q 2 . Fig.Q shows the definition of the particle momenta and helicities. The 
momentum transfer (squared) is given by A = P' - P (t = A 2 ). Using the light 
cone vector n we define the longitudinal momentum transfer by £ = —n ■ A/2. 
The proper definition of the twist- 2 tensor or quark helicity flip GPDs Ht, Et, 
Ht and Et can be found in 4. GPDs provide a solid framework in QCD to 
relate many different aspects of hadron pyhsics: 

• The forward limit A — > of certain GPDs reproduces the well known 
parton distributions, that is H(x, £ = 0, t = 0) = q(x) = fi{x), H(x, 0, 0) = 
Aq(x) = g\(x) and Ht(x,0, 0) = 5q(x) = h\(x). 

• The integral over the longitudinal momentum fraction j dx of the GPDs 
gives the Dirac, Pauli, axial, pseudo-scalar, tensor etc. form factors, 
/ dxH{x,£,t) = Fi(t), f dxE = F 2 (t), f dxH = g A (t), f dxE = g P (t), 
J dxHx = gr(t) etc. 

• The Fourier transforms / dA ± e~ ib± - A± of the GPDs H, H and H T at £ = 
are coordinate space probability densities in the impact parameter b± [JJ]. 

• The forward limit of the x-moment of the GPD E, J dxxE(x, 0, 0) = 
B20 (0), allows for the determination of the quark orbital angular momen- 
tum contribution to the nucleon spin, L q = l/2((x) + -B20 — Aq), where (x) 
is the quark momentum fraction. 

The fact that the impact parameter dependent quark distributions like 



have the interpretation of probability densities in the transverse plane marks 
a major step forward in the understanding of the nucleon structure. A very 
interesting observation in this respect is that at large longitudinal momentum 
fraction, x — > 1, i.e. when the active parton carries the whole nucleon momentum 
and the spectators give only a negligible contribution, the active parton becomes 
the center of momentum of the nucleon at R± = 0. Accordingly, the impact 
parameter profile of the nucleon will be strongly peaked at b± = 0. Keeping 
in mind that the parton distributions eventually vanish in the limit x — ► 1, we 
expect that lim^i q(x, b±)/q(x) oc 5 2 (b±). Transferred to momentum space, we 
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Figure 2. Some parameters of our calculation. 



have the striking prediction that the normalized GPD H(x,0,t)/q(x) will be 
constant in the limit x — > 1 [HI El • 

Interesting recent results on GPDs using available experimental data on form 
factors and PDFs can be found in [Hj. 

On the lattice side, we cannot deal directly with matrix elements of bilocal 
light-cone operators. Therefore we first transform the LHS in Eqs. (|H2[) to Mellin 
space by integrating over x, i.e. J dxx n ~ l . This results in matrix elements of 
towers of local operators which are in turn parameterized in terms of so-called 
generalized form factors (GFFs) A(t), B(t) . . ., 

(P',\'\$(0)riD -iD V(0)|P,A) = U(P',X')(a l p- t " n A(t) 

+ b^-^B{t) + ---)U{P,\),{A) 

where the subtraction of traces is implicit. The explicit parameterization for 
the tower of vector operators, r = 7^, is given in [S| in terms of the GFFs 
A n i(t) , B n i(t) and C n o(t). In the axial vector case, the corresponding expression 
in terms of the GFFs A n i (t) and B n i (t) is shown in ^01 ■ Parameterizations for the 
tensor GPDs have been derived recently in [111 112] . Having computed the GFFs 
in e.g. lattice QCD, it is fairly simple to get the moments of the corresponding 
GPDs by using the following polynomial relations 13 

1 n—l 
H n {^t) = J ' dxx n - 1 H{x^t)= £ (-2O l A ni (t)+(-2O n C n0 (t)\ neven , 

_^ i=0,even 
n-1 

E n {^t) = (-^yB ni (t)-(-20 n C n0 (t)\ neven , 

i=0,even 

n—l n—l 

H n (Z,t) = Yl i-XYAuV), E n (£,t)=Y (-^YBniit). (5) 

i=0,even i=0,even 

In the limit £ — ► 0, the complete nucleon structure related to H and H is accord- 
ing to Eq.fJSJ) represented by the sets of GFFs A n o(t) and A n o(t), n = 1 ... 00. 



2 Lattice results 

The lattice calculation and extraction of the GFFs from nucleon two- and three- 
point-functions follows the standard procedure as outlined in e.g. [TH ITdI ITB] . 
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Figure 3. Lowest three moments of the GPDs H(x, 0, t), E(x, 0, i) (left column) and H(x, 0, t), 
E(x,Q,i) (right column). 



Our dynamical computations are currently based on a set of coupling- and quark 
mass-combinations with up to 0(1200) configurations and with pion masses rang- 
ing from 1.1 down to 0.55 GeV. The lattice spacing lies in between ~ 0.07 
and ~ 0.12 fm. Since we are using non-perturbatively 0(a) improved Wilson 
fermions, we do not expect strong discretization effects |16j . A list of the param- 
eters is given in Fig.(j2j). We do not take into account the disconnected contribu- 
tions which are numerically much harder to evaluate. This is irrelevant as long as 
we concentrate on the iso- vector channel u — d (as we do here), for which the dis- 
connected pieces would cancel due to iso-spin symmetry. To give an overview of 
our results we show in Fig.© the lowest three moments of the GPDs H(x,0,t), 
E(x,0,t) and H(x,0,t) as well as the lowest two moments of E(x,0, t). Up to 
now, only the lowest moments n = 1 which we show here are non-perturbatively 
renormalized, and operator mixing has to be taken into account for a proper 
renormalization of the x 2 -moments for i > |17| . We do not expect the renor- 
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Figure 4. Change of slope of the GFFs (left)and narrowing of the corresponding b± dependent 
GFFs (right); n = 1 corresponds to the lowest (widest) curve in the left(right)-hand plot. 



malization to affect the following general discussion and the extraction of the 
dipole masses below. Please note that for reasons of clarity i?3o(t) has been 
shifted upwards by 0.05, as indicated by the superscript arrow in Fig.©, lower 
left panel. All GFFs have been fitted using a standard dipole ansatz, 

<r le (t) = , (6) 

(1 - t/m 2 D ) 

except for the pseudo-scalar FF B w (t) = gp(t), where we used a dipole ansatz 
including the pion-pole, 

i?r o pole (t) = 10 2 [ ' . (7) 

{l-t/mlf(l-t/ml) 

The fits given by the gray bands in Figs.® show an overall satisfactory de- 
scription of the lattice numbers. While the statistical errors for the lowest two 
moments are rather small, they tend to be somewhat larger for the x 2 -moment. 
The statistics are particularly poor for ^32, ^32 and the B%i GFFs, which are not 
shown in this work. It turns out that the moments of H and H in the iso- vector 
channel are a factor of 2 ... 5 smaller than the corresponding moments of E and 
E, while C20 which contributes to both H and E is compatible with zero. Our 
results indicate furthermore strong cancellations of u and d flavor contributions 
in the iso-singlet channel for the helicity independent GPD E u+d . 

Let us now take a look at the i-dependence of the GPDs and the impact 
parameter dependent quark distribution q(x,b±). The observation from section 
(^Q) that lim x ^i H(x, 0, t)/q{x) = 1 translates directly to the statement that the 
GFFs A n Q(t) become constant for very large n since increasingly more weight 
lies on values of x close to one for lim^^oo. To be precise, we expect that 
linin^oo A n o(t) / A n o(0) = 1, independent of t. Since the values for the momen- 
tum transfer squared in our computation are unevenly distributed in the interval 
—3.5 < t < GeV, and since there is in particular a gap between t = and the 
smallest non-zero t which can be realized on the given lattices, t m i n , it is not 
useful to perform a discrete Fourier transformation (FT) of our results to impact 
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Figure 5. Joined results for all /3, n combinations for the lowest two moments of H (x, 0, t) (left 
column) and H(x,0,t) (right column). 



parameter space. Instead, we directly use the continuous FT of the dipole ansatz 
Eq.©. The normalized impact parameter dependent moments q n=1 ' 3 (b±) with 
J db±q n=1 '" S (b±) = 1 are shown in Fig.(j3J. For reasons of comparison we show 
in Fig.(jlJ) also the corresponding GFFs, normalized to ^^^(O) = 1. The gray 
bands in the impact parameter plot reflect the errors associated with the dipole 
masses, Amo- The results in Fig.Q show very nicely the anticipated flattening 
of the GFFs and the narrowing of the corresponding distributions in impact pa- 
rameter space, going from n = 1 to n = 3. This flattening of the GPDs has first 
been observed in [T5] . 

Since we have results for a large number of pion masses at hand, it is manda- 
tory to concentrate in the following on the important issue of the dependence 
and the chiral extrapolation of our calculation. Due to the use of Wilson fermions, 
the dominant part of our computation is for pion masses above 600 MeV where 
chiral perturbation theory in its current form is most probably not applicable. In 
particular for the dipole masses mp, we therefore resort to a purely phenomeno- 
logically motivated pion mass dependence of the form = m® D + arn^.. We 
have checked that a linear dependence on m n is clearly unfavored by our re- 
sults. For the following analysis we have normalized our results to A^ olc (0) = 1 
first. One option is then to perform dipole fits for all individual /3, k combi- 
nations for a given moment n and fit the pion mass dependence of the dipole 
masses with the given ansatz. However to optimally utilize the information on 
the dipole masses from our collection of results, it is better to use a combined 
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depending on the two parameters mP D and a. Having done the fit, we may shift 
the raw numbers to a common curve given by Eq.@ at the physical pion mass, 
m 7r,phys = 139 MeV. The results of this procedure are shown in Fig.© for the 
GFFs A n o and A n Q, n = 1,2. The fit parameters are found to be 
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Although the corresponding rather large error does not allow for a definite con- 
clusion, it is interesting to note that the slope a changes sign going from A\q to 
A20. We observe an overall nice clustering of the lattice points in Fig.© much 
in favor of our ansatz Eq.(jHJ). The actual fits are given by the gray bands and 
show a very good x 2 /NODOF lying below 0.2 for the helicity independent and 
below 0.5 for the helicity dependent case. In particular for A n Q, a small fraction 
of the lattice points seem to be significantly below the main line. Future investi- 
gations have to show if this is due to an unaccounted dependence on the lattice 
spacing, the (finite) lattice volume or just insufficient statistics. An attempt to 
compare our results for the axial and tensor coupling, 5^4(0) and st(0), with 
chiral perturbation theory is described in ^H]- We plan to extract the quark or- 
bital angular momentum contribution to the nucleon spin as soon as the missing 
renormalization constants become available. 



3 Conclusion 

We have computed the lowest moments of GPDs for a number of different pion 
masses in the range above 550 MeV. The results are very well described by a com- 
bined dipole-mass/pion-mass fit. While we observe no sign of chiral logarithms 
in the available range of pion masses, let us note that first promising lattice cal- 
culations towards the chiral region are currently under way using overlap t 20 as 
well as domain wall fermions |21j . 
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